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RESUME 


Nous  decrivons  une  technique  experimental  permettant  de  deter¬ 
miner  le  seuil  d ' endommagement  des  surfaces  optiques  de  verre  et  de 
plastique  irradiees  par  un  laser  CC^  TEA  impulsionnel .  Les  densites 

d'energie  correspondant  5  une  probability  d ' endommagement  de  50%,  telles 

2 

qu'estimees  par  une  analyse  probit,  sont  de  4.0  J/cm  pour  la  si lice 

2  2 
fondue,  de  2.1  J/cm  pour  le  verre  sodocalcique ,  de  2.3  J/cm  pour  le 

2 

thermoplastique,  de  2.3  J/cm  pour  le  plastique  acrylique  et  de 
2 

1.7  J/cm  pour  le  polystyrene.  Nous  presentons  une  interpretation  des 
distributions  statistiques  observees  en  fonction  de  divers  mecanismes 
d ' endommagement .  Enfin,  plusieurs  ameliorations  susceptibles  d'accroT- 
tre  la  precision  des  resultats  obtenus  sont  suggerees.  (NC) 


ABSTRACT 

An  experimental  technique  is  described  for  determining  the 
damage  thresholds  of  glass  and  plastic  optical  material  surfaces, 
opaque  to  laser  radiations,  when  irradiated  with  a  pulsed  TEA-Cd^ 
laser.  Use  of  a  probit  analysis  to  estimate  the  irradiance  required 

for  a  50%  probability  of  damage  gives  threshold  levels  of  4.0  J/cm^  ' A  ’  S 

2  2 
for  fused  silica,  2.1  J/cm  for  soda-lime  glass,  2.3  J/cm  for 

2  2 

polycarbonate,  2.3  J/cm  for  acrylic  and  1.7  J/cm  for  polystyrene 
plastics.  The  significance  of  the  statistical  distributions  observed 
in  attributing  various  mechanisms  to  the  damage  is  discussed.  Finally, 
some  improvements  are  suggested  to  make  the  measurements  more  accurate .' (U) 
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1.0  INTRODUCTION 

A  number  of  workers  have  shown  that  damage  to  transparent  optical 
materials  from  laser  radiation  is  a  statistical  phenomenon.  Bass  and 
Barret  (Ref.  1)  first  demonstrated  the  statistical  nature  of  the  damage 
due  to  an  intrinsic  electron  avalanche  breakdown,  while  others  (Refs. 
2-5)  showed  that  damage  due  to  absorbing  inclusions  and  structural 
defects  in  the  material  were  also  statistical.  In  these  cases,  the 
concept  of  a  damage  threshold  could  only  be  defined  statistically  since 
irradiating  apparently  identical  sites  with  identical  laser  pulses 
produced  damage  for  some  of  the  trials  only. 

On  the  other  hand,  homogeneous  absorption  by  the  irradiated 
material  is  non  statistical,  and  materials  damaged  by  this  phenomenon 
should  show  a  definite  threshold.  Under  this  assumption,  the  absence 
or  presence  of  statistics  could  be  an  indication  of  the  damaging  mech¬ 
anism.  Caution  is  required,  however,  since  there  are  many  sources  of 
statistical  fluctuations  associated  to  damage  threshold  experiments. 

For  example,  lasers'  inherent  pulse-to-pulse  variation  in  energy  may 
mask  the  statistics  of  the  damage  phenomenon;  the  surface  conditions 
of  the  samples  may  influence  the  amount  of  energy  absorbed;  the  actual 
failure  mechanism,  triggered  by  the  absorbed  radiation,  may  be  statisti¬ 
cal  even  if  the  absorption  is  not. 

In  this  work,  we  present  some  typical  results  of  damage  experi¬ 
ments  on  various  glasses  and  plastics,  opaque  to  the  laser  radiation. 
Damage  is  defined  on  the  basis  of  changes  which  occur  in  the  physical 
appearance  of  the  surface  of  the  materials.  Since  the  absorption  is 
homogeneous  in  the  materials  chosen,  the  damage  is  presumed  to  be  non- 
statistical . 
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To  measure  distributions  in  the  damage  thresholds,  we  adapted 
the  experimental  technique  used  by  Hacker  and  Halverson  (Ref.  6)  for 
determining  the  breakdown  threshold  of  gases.  The  influence  of  laser 
energy  fluctuations  is  largely  eliminated  by  recording  the  energy  of 
each  pulse  along  with  an  indication  of  whether  or  not  damage  occurred. 
Any  remaining  statistics  that  can  be  attributed  to  the  damage  phenom¬ 
enon  may  be  used  to  characterize  its  process. 

Section  2  describes  the  experimental  setup  and  the  procedure  for 
acquiring  the  data.  An  outline  of  the  statistical  analysis  which  allows 
to  estimate  the  damage  threshold  and  the  parameters  of  the  probability 
density  function  appears  in  Section  3.  The  results  and  a  discussion  of 
their  significance  follow  in  Section  4. 

This  work  was  performed  at  DREV,  between  January  1977  and  May 
1980,  as  partial  fulfilment  of  tasking  for  DLAEEM  under  PCN  33B39, 

Laser  Induced  Damage  to  Optical  Surfaces. 

2.0  EXPERIMENTAL 
2 . 1  Damage  Testing  Facility 

Figure  1  schematically  illustrates  the  physical  setup  used  for 
the  damage  tests.  Damage  was  induced  by  a  Laflamme  type  (Ref.  7) 

TEA-CO,  laser  (Ref.  8)  modified  for  variable  pulse  length  operation 
(Ref.  9).  The  pulse  length  of  this  laser  can  be  varied  from  0.1  us 
to  100  us  (width  at  10%  maximum  height)  by  changing  the  reflectivity 
of  the  output  coupler  and  the  composition  of  the  gas  mixture.  Figure  2 
shows  the  typical  waveform  of  the  34 -us  pulses  used  in  these  experi¬ 
ments.  The  available  pulse  energy  ranges  from  5  to  50  J,  depending  on 
the  pulse  length  and  the  discharge  voltage.  The  laser  operation  is 
multimode  on  the  P20  line,  at  10.6  um,  with  a  beam  approximately  50  mm 
by  50  mm  at  the  output. 
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FIGURE  1  -  Physical  setup  for  laser  damage  threshold  experiments.  The 
pressure  of  propylene  in  the  gas  cell  attenuator  controls 
the  beam  energy.  The  two  reflected  beams  from  the  NaCl,  5°- 
wedge  beamsplitter  are  monitored  by  the  pyroelectric  Joule- 
meter  for  a  measure  of  the  beam  energy,  and  by  the  photon 
drag  detector  for  the  time  evolution  of  the  laser  pulse.  The 
spherical  mirror  and  the  optical  integrator  produce  a ~ 12  mm 
square  spot  on  the  target  plane.  Long-focal -length  mirrors 
were  used  to  avoid  the  possibility  of  air  breakdowns  in  the 
focal  regions.  Except  at  the  flat  turning  mirror,  all  angles 
between  incident  and  reflected  beams  were  kept  less  than  10° 
to  reduce  polarization  effects  at  the  surface  of  dielectrics, 
and  aberrations  from  the  spherical  mirrors. 
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FIGURE  2 
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-  Time  evolution  of  the  laser  pulse  recorded  from  a  photon  drag 
detector  with  a  150 -MHz -bandwidth  oscilloscope  terminated 
with  SO  fi.  Pulse  width  of  34  ys  is  measured  as  the  width  at 
10%  of  the  maximum. 
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A  series  of  plane  and  spherical  mirrors  (Fig.  1)  folded  and 
focussed  the  beam  to  a  spot  approximately  12 -mm  square  at  the  target 
plane.  A  homogenizing  reflector,  consisting  of  32  small  flat  mirrors 
mounted  on  a  spherical  surface,  was  also  incorporated  to  obtain  a 
uniformly  distributed  beam  profile.  This  arrangement  avoided  air 
breakdown,  and  allowed  the  use  of  small  off-axis  angles  to  reduce 
the  effects  of  spherical  aberrations.  It  also  provided  for  spatial 
filtering  if  a  more  uniform  beam  profile  were  needed  for  greater  pre¬ 
cision  in  future  experiments.  The  current  profile  at  the  target  plane 
is  flat  to  within  10%.  Provision  was  made  for  a  second  NaCl,  5°-wedge 
beam  splitter  to  be  inserted  between  the  final  focussing  mirror  and 
the  target  plane  to  provide  two  additional  sample  beams,  one  for 
monitoring  the  laser  wavelength  by  means  of  a  C 0^  laser  spectrum  ana¬ 
lyzer,  the  other  for  future  continuous  monitoring  of  the  laser  beam 
profile  at  the  target  plane  with  a  pyroelectric  vidicon  camera. 

A  mechanism  (Fig.  3)  held  the  target  and  repositioned  it  between 
the  shots  to  expose  a  fresh  site  to  the  laser  beam.  This  mechanism 
allowed  the  laser  beam  to  scan  the  target  in  raster  fashion  to  cover 
rectangular  samples  from  5  cm  to  30  cm  on  a  side. 

2 . 2  Procedure 


The  technique  used  was  a  multiple-shot,  1-on-l  experiment  in 
which  each  site  of  the  sample  was  irradiated  once  by  a  pulse  from  the 
CO 2  laser.  This  avoided  cumulative  effects,  such  as  laser  polishing 
or  conditioning  of  the  sample  (Ref. 10),  and  ensured  that  the  effect  of 
each  shot  was  independent.  To  obtain  as  many  damaging  shots  as 
undamaging  ones,  we  controlled  the  laser  energy  from  shot-to-shot 
with  a  gas  cell  attenuator.  We  attempted  to  collect  sufficient  data 
at  each  energy  level  to  ensure  at  least  one  surviving  or  one  damaged 
site . 


UNCLASSIFIED 

(> 


CURE  3  -  Mechanical  scanner  used  to  transport  a  fresh  sample  site 
into  position  between  laser  shots.  A  5  to  50  cm  square 
sample  may  be  accommodated;  it  moves  horizontally  between 
adjustable  limits  then  shifts  20  mm  vertically  to  start  a 
new  horizontal  scan.  Scan  speed  may  be  set  over  a  wide 
range  hut,  in  this  experiment  it  was  adjusted  to  space  sample 
sites  ~20  mm  apart  between  laser  shots. 
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The  joulemeter  produces  a  peak  voltage  pulse  proportional  to 
the  laser  beam  energy  it  monitors.  The  corresponding  energy  actually 
falling  on  the  target  plane  can  be  calculated  if  the  ratio  between  the 
joulemeter  voltage  and  the  target  energy  is  known.  This  ratio  was 
obtained  by  direct  measurement,  and  frequently  verified  so  as  to  detect 
any  variation  caused  by  inadvertent  changes  in  the  alignment  or  the 
optical  properties  of  the  components. 

To  increase  precision,  we  offset  the  joulemeter  voltage  by  a 
known  constant  before  reading  its  peak  with  a  digitizer  and  storing 
it  in  one  of  two  256-channel  memories.  Short  duration  noise,  which 
appears  on  the  joulemeter  cable  at  the  time  of  the  laser  firing, 
prevents  the  use  of  the  peak  reading  mode  of  the  digitizer.  However, 
since  the  joulemeter  voltage  rises  slowly,  over  several  milliseconds, 
and  reaches  its  maximum  after  the  noise  has  subsided,  the  reading 
could  be  taken  by  triggering  the  digitizer  with  a  pulse  delayed  to 
coincide  with  the  peak. 

Damage  was  detected  visually  since  all  the  sample  materials  were 
transparent  in  the  visible  region  of  the  spectrum.  The  polystyrene, 
polycarbonate  and  silica  surfaces  were  damaged  by  melting.  Damage  in 
these  materials  was  detected  by  projecting  images  of  their  surfaces  on 
a  screen,  with  a  white  light,  to  observe  their  patterns.  It  appeared 
as  white  spots  in  the  acrylic  plastic  (PMMA) ,  and  as  fine  cracks,  which 
occurred  within  5-10  s  of  irradiation,  in  the  glass.  In  these  two 
materials,  damage  was  detected  by  illuminating  the  samples  with  a 
collimated  beam  of  white  light  from  an  arc  lamp  and  by  viewing  the 
laser  irradiated  area  with  a  standard  TV  monitor. 


»  , 
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A  555  timer,  connected  as  a  bistable  multivibrator,  registered 
whether  a  site  was  damaged  or  not.  When  damage  occurred,  the  multi¬ 
vibrator  output  was  manually  set  high,  otherwise,  it  was  left  low. 

The  pulse  that  triggered  the  digitizer  for  the  joulemeter  also  trig¬ 
gered  a  second  digitizer  which  recorded  the  state  of  the  multivibrator 
in  the  second  256-channel  memory  and  reset  it  low.  This  allowed 

several  seconds,  after  the  laser  shot,  for  the  damage  decision  to  be 

t  h 

made.  We  must  be  aware  that,  although  the  energy  of  the  n  laser 

t  Vi 

shot  was  recorded  in  the  n  channel  of  the  first  memory,  the  damage 

t  h 

decision  was  recorded  later  in  t,  the  n+1  channel  of  the  second 
memory.  This  can  easily  be  integrated  to  the  computer  program  to 
perform  the  analysis. 

Via  a  CRT  terminal,  the  data  computed  from  the  numerous  shots 
was  transferred  from  the  memories  to  DREV's  main  computer  through  a 
custom-built  interface.  Within  a  few  seconds,  preliminary  results  of 
the  statistical  analysis  were  displayed  on  the  CRT  terminal's  screen 
(Figs.  4-8)  so  that  one  could  rapidly  determine  if  the  data  was 
adequate . 

The  samples  were  tested  as  received  from  the  distributors.  No 

special  treatments  were  applied  to  the  surfaces  beforehand  other  than 

the  removal  of  any  obvious  dirt  and  dust  with  compressed  air.  The 

soda-lime  glass  was  ordinary  window  glass  from  the  carpentry  shop. 

R* 

The  polycarbonate  plastic  sheet  was  Lexan  9030,  the  acrylic  plastic 
TM** 

was  Plexiglas  ,  and  the  polystyrene  was  a  "K-LUX"  Safe-t-vue  clear 
sheet.  The  silica  was  a  Corning  7940  fused-silica  plate  with  a 
commercial  polish. 

*  Trademark  of  the  General  Electric  Company,  Pittsfield,  Mass.  01201 
TM 

**  Trademark  of  the  Rohm  and  Hass  Co.,  Philadelphia  Pa.  19105 
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FIGURE  4  -  Experimental  results  for  the  measurement  of  the  damage 

threshold  of  fused  silica  irradiated  by  a  34-ys-long  TEA-CO, 
laser  pulse  at  10.6  um.  The  upper  left  graph  is  a  histogram 
of  all  laser  shots  as  a  function  of  laser  energy  (represented 
by  the  signal  voltage  from  the  energy  meter).  The  lower  left 
graph  is  a  histogram  of  those  shots  which  damaged  the  sur¬ 
face.  (In  fused  silica,  the  damage  was  characterized  by  the 
onset  of  surface  melting).  The  upper  right  graph  shows  the 
probability  of  causing  damage  as  a  function  of  laser  energy. 
The  crosses  represent  the  relative  frequencies  of  damage 
obtained  by  dividing  the  histogram  of  damaging  shots  by  that 
of  the  total  number  of  laser  shots.  The  curve  is  the  max¬ 
imum  likelihood  estimate  of  the  probability  where  a  normal 
probability  density  distribution  was  assumed.  The  lower 
right  hand  graph  is  a  probit  plot  of  the  same  data.  The 
probit  is  a  coordinate  transformation  which  linearizes  the 
probability  curve.  The  dashed  lines  represent  95%  confidence 
limits  based  on  Student's  t  distribution.  The  damage  thresh¬ 
old  corresponding  to  the  50%  probability  of  damage  occurs  at 
0.120  t  0.002  V,  which  represents  an  irradiancc  of  4.0  .I/cm". 
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FIGURE  5  -  Experimental  results  for  the  measurement  of  the  damage 

threshold  in  soda- lime  glass.  An  explanation  of  the  graphs 
appears  in  the  caption  of  Fig.  4.  In  soda-lime  glass,  damage 
is  characterized  by  surface  cracking.  The  threshold  for  a 
50%  probability  of  damage  is  0.242  ±  0.003  V,  corresponding 
to  an  irradiance  level  of  2.1  J/cnT. 
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Signal  CVolfc-s)  Signal  CVo I ts) 

FIGURE  6  -  Experimental  results  for  the  measurement  of  the  damage 

threshold  in  polycarbonate  plastic.  An  explanation  of  the 
graphs  appears  in  the  caption  of  Fig.  4.  In  polycarbonate 
damage  is  characterized  by  surface  melting.  The  threshold 

for  a  50%  probability  of  damage  is  0.075  ±  0.001  V,  corre- 

? 

sponding  to  an  irradiance  level  of  2.3  J/cm  . 


Damaged  s i tes  Laser  shots 
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N  3  POLYSTYRENE  PLASTIC  VED  16  APR  1980 


Signal  CVoltsJ  Signal  CVollsJ 

FIGURE  8  -  Experimental  results  for  the  measurement  of  the  damage 


threshold  in  polystyrene.  An  explanation  of  the  graphs 

appears  in  the  caption  of  Fig.  4.  In  polystyrene,  damage 

is  characterized  by  surface  melting.  The  threshold  for  a 

50%  probability  of  damage  is  0.0579  ±  0.0009  V,  corresponding 

2 

to  an  irradiance  level  of  1.7  J/cm  . 
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3.0  STATISTICAL  ANALYSIS  OF  DAMAGE.  MEASUREMENTS 


3.1  Sources  of  Statistics  in  Threshold  Measurements 


The  five  materials  chosen  for  this  series  of  experiments  are 
strong  absorbers  of  the  10.6-ym  radiation  from  CO.,  lasers.  Assuming 
that  the  absorption  is  homogeneous,  one  at  first  expects  the  damage 
due  to  the  absorbed  energy  to  be  deterministic,  i.e.  free  from  statis¬ 
tics.  A  few  measurements  quickly  showed  that  this  is  not  strictly 
true  as  two,  apparently  similar,  laser  shots  produced  different 
results.  The  statistical  nature  of  these  results  could  arise  from 
several  causes  such  as  fluctuations  in  the  laser  output,  lack  of 
material  homogeneity,  presence  of  localized  dirt  or  polishing  compounds, 
variations  in  the  positioning  of  the  sample,  fluctuations  of  the  laser 
beam  profile,  or  the  intrinsic  way  in  which  damage  manifests  itself 
after  the  laser  energy  is  absorbed. 

This  experiment  was  designed  to  eliminate  the  effects  of  the 
variations  in  the  laser  output  by  accurately  measuring  the  energy  of 
each  pulse.  Lack  of  material  homogeneity  seems  unlikely  in  our 
experiments  since  the  sites  were  all  on  the  same  sample,  and  any  small 
variations  in  the  absorption  are  negligible  compared  to  the  90%  to 
96%  absorption  (Ref.  12)  for  these  materials.  Although  the  samples 
were  tested  as  received,  their  appearance  ruled  out  any  gross  contami¬ 
nation  by  foreign  substances.  In  any  case,  the  absorption  is  not 
expected  to  be  greatly  affected  by  small  amounts  of  imperceptible 
dirt.  Variations  in  the  positioning  of  the  sample  were  reduced  by 
ensuring  that  the  scanning  motion  was  parallel  to  the  surface  of  the 
sample.  A  poor  alignment  of  the  scanner  would  have  caused  the  damaged 
sites  to  collect  in  the  same  portion  of  the  sample.  Temporal  variations 
in  the  laser  pulses  were  not  expected  to  affect  results  since  those 
pulses  were  all  nominally  the  same  length,  and  any  small  variations 
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occurred  over  time  periods  much  shorter  than  the  diffusion  time 
necessary  for  the  irradiated  areas  to  cool  down.  However,  fluctuations 
in  the  beam  profile  and  in  the  way  damage  occurs  in  glass  do  affect 
the  results,  as  discussed  in  Section  4.  Once  all  statistics  due  to 
extrinsic  causes  have  been  accounted  for,  those  that  remain  must  be 
attributed  to  intrinsic  ones  and  they  must  be  checked  against  their 
exp ec ted  d i s t ribut i on . 

3 . 2  Probit  Analysis 


The  underlying  assumption  is  that  the  probability  for  damage 
increases  with  the  density  of  energy  falling  on  the  surface  of  the 
material.  Furthermore,  we  assume  that  the  threshold  levels  are 
normally  distributed,  with  a  variance  related  to  the  damage  mechanism. 
If  the  energy  density,  or  irradiation,  is  denoted  by  E,  the  probability 
distribution  of  damaged  sites  may  be  expressed  by 

dP  =  (27ra2)"1/2  exp  [- (E-y)2/2o2]dE  [1] 

where  dP  is  the  probability  of  damage  occuring  over  the  energy  density 

2 

range  dE,  y  is  the  mean  of  the  distribution,  and  a  is  the  variance. 

The  analysis  used  to  estimate  the  parameters  u  and  of  eq.  1 
is  based  upon  the  probit  transformation  of  the  experimental  results 
(Ref,  13) .  The  probit  is  the  random  variable  which  corresponds  to  the 
probability  in  a  normal  distribution,  with  mean  5  and  variance  1. 
Symbolically,  the  probit,  Y,  of  the  probability,  P,  is  defined  by 
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P  =  (2tt) 


-1/2 


/'S 


exp  [-u  /2]  du 


[2] 


Integration  of  eq.  1  gives  the  expected  proportion  of  samples  damaged 
up  to  an  energy  density  Eq  as  ‘ 


P  = 


(2raV1/2 


exp  [-(E-y)2/2o2]  dE 


[3] 


Comparison  of  [2]  and  [3]  shows  that  the  probit  of  the  expected 
damaged  proportion  is  related  to  the  energy  density  by  the  linear 
equation 


Y  =  5  +  (Eq  -  y)/a  [4] 

Traditionally,  the  probit  method  is  to  fit  a  best  straight  line 
to  eq.  4,  by  adjusting  y  and  a  to  obtain  the  maximum  likelihood 
estimate  for  the  experiment.  That  is,  the  two  adjustable  parameters 
of  the  distribution  function  are  varied  to  obtain  the  probability  dis¬ 
tribution  which  corresponds  to  the  results  observed  to  be  the  most 
likely  outcome  of  the  experiment.  In  the  past,  the  probit  transfor¬ 
mation  was  used  to  linearize  the  probability  distribution  and  thus 
to  facilitate  this  fitting.  Although  modem  computers  makes  this  no 
longer  strictly  necessary,  expressing  the  results  in  terms  of  the 
probit  is  useful  for  comparison  and  interpretation.  Our  fitting 
procedure,  which  closely  follows  the  method  developed  in  Appendix  II 
of  Ref.  13,  is  performed  in  APL,  as  described  in  Appendix  A.  Once 
the  maximum  likelihood  estimate  of  the  probability  distribution 
function  has  been  obtained,  it  is  easy  to  determine  the  energy  density 
required  for  a  given  probability  of  damage.  For  a  50%  probability  of 
damage,  the  energy  level  required  is  equal  to  y. 
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Strictly,  the  probit  refers  only  to  the  normal  distribution. 
However,  it  is  instructive  to  perform  a  similar  analysis  for  a  general 
distribution.  In  particular,  the  APL  program  we  elaborated  can  also 
perform  the  analysis  for  a  uniform  distribution. 

4.0  RESULTS  AND  DISCUSSION 

With  an  APL  program,  we  scaled  and  sorted  the  raw  data  into 
two  histograms,  one  of  the  total  number  of  laser  shots,  the  other  for 
those  shots  causing  damage  as  a  function  of  joulemeter  voltage.  Fol¬ 
lowing  the  probit  analysis,  the  data  was  displayed  graphically,  as 
shown  in  Figs.  4  to  8,  for  fused  silica,  soda-lime  glass,  polycarbonate, 
acrylic  and  polystyrene  plastics  respectively.  In  each  of  those 
figures  the  upper  left  graph  is  the  histogram  of  all  laser  shots 
plotted  as  a  function  of  the  joulemeter  signal  voltage  proportional 
to  the  energy  density  of  the  beam  at  the  target.  The  lower  left  graph 
is  the  histogram  of  those  shots  which  damaged  the  surface.  The  upper 
right  graph  shows  the  probability  of  causing  damage  as  a  function  of 
the  joulemeter  voltage.  The  crosses  represent  the  relative  frequencies 
of  damage  obtained  by  dividing  the  histogram  of  damaging  shots  by  that 
of  the  total  number  of  laser  shots.  The  curve  is  the  maximum  like¬ 
lihood  estimate  of  a  normal  probability  density  distribution.  The 
lower  right  graph  is  a  probit  plot  of  the  same  data.  The  dashed  lines 
represent  the  95%  confidence  limits  based  on  Student's  t  distribution. 

We  define  the  damage  thresholds  as  the  energy  density  levels 
at  which  the  probability  of  damage  is  50%.  Table  I  gives  the  values 
obtained  for  the  thresholds,  and  the  standard  deviation  of  the  proba¬ 
bility  distribution  function. 
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MAHKIAI. 

THRhSHOLl)  (uj 

STANDARD  III  V!  ATI  ON  (■->) 

COU  FIC  HAT  OF 

VARIATION 

l  .l/cm" ) 

fd/cm" 

i 

Fused  silica 

4.0  (.05)* 

0.31 

( 0 . 0b 1 

0.076 

f  .016) 

Soda- lime  glass 

J.l  (.01) 

0.089 

(0.01) 

0.042 

f .  air ) 

Polycarbonate  plastic 

2.3  (.02) 

0.16 

(0.03) 

0.071 

(  .013) 

Acrylic  plastic 

2.3  (.01) 

0.047 

(0.01) 

0 . 02 1 

r.oo5i 

Polystyrene  plastic 

1.’  (.01) 

0.091 

(0.01) 

O.OS4 

( .0081 

*  Numbers  in  parentheses  are  the  standard  deviations  of  the  estimates  for  the  parameters 
of  the  probability  distribution  function. 
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FIGURE  9  -  Analysis  of  the  experimental  results  for  acrylic  plastic 
assuming  a  uniform  density  for  the  probability  density 
function.  A  full  explanation  of  the  graphs  appears  with 
Fig.  4.  The  estimated  energy  density  for  a  50%  probability 

of  damage  occurs  at  0.304  ±  0.002  V,  corresponding  to  an 

2 

energy  density  of  2.3  J/cm  .  This  result  differs  insignif¬ 
icantly  from  that  obtained  using  a  normal  distribution 
function . 
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Figure  9  illustrates  the  results  of  the  soda-lime  glass  analysis 
when  a  uniform  probability  distribution  is  used.  The  results  for  a 
and  y  only  slighty  differ  from  those  of  the  normal  distribution  in 
Fig.  5.  There  is  no  significant  difference  between  the  two  distribu¬ 
tions  for  representing  the  experimental  data.  This  holds  for  the  other 
materials  as  well . 

The  coefficient  of  variation,  (a/u) ,  is  an  estimate  of  the 
relative  dispersion  in  the  threshold  data.  As  shown  in  Table  I,  this 
coefficient  is  small  for  all  materials  tested.  The  finite  binwidths 
used  in  calculating  the  estimates  had  no  significant  effect  since 
Sheppard's  correction  was  negligible  and  recalculation  of  the  estimates 
of  a  showed  no  changes,  even  when  bins  several  times  larger  were  used. 

The  experiments  on  the  fused  silica,  polycarbonate,  and  poly¬ 
styrene  samples  yielded  the  largest  coefficients  of  variation.  Since 
it  meant  looking  for  surface  distortions  due  to  melting,  detecting 
damage  in  these  three  materials  was  difficult.  In  fact,  their  melting 
points  are  ill-defined  and  the  recorded  deflection  and  heat  resis¬ 
tance  temperatures  factors  for  plastics  vary  greatly  (Ref.  14) .  Some 
of  the  variations  in  the  reported  temperatures  can  be  attributed  to 
differences  between  manufacturing  processes  or  batches.  On  the  other 
hand,  silica  is  a  true  glass;  it  does  not  'melt',  but  because  of  a 
decreasing  viscosity,  it  becomes  softer  with  increasing  temperature. 

However,  the  cracks  in  the  glass  and  the  white  spots  in  the  PMMA 
samples  were  readily  seen  and  gave  a  clear  indication  of  the  damage 
occurence.  Therefore,  the  dispersions  observed  in  the  fused  silica, 
polycarbonate  and  polystyrene  seem  largely  due  to  the  damage  detecting 
method.  The  observations  in  PMMA  and  glass  likely  reflect  other 


causes . 
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The  damage  in  PMMA  consisted  of  many  voids,  ~10  urn  in  size, 
lying  in  the  surface  (Fig.  10).  There  is  little  evidence  of  ejected 
debris.  This,  and  the  strong  solvent-like  odor  exhaled  during  the 
experiments  suggest  that  damage  occured  when  the  plastic  decomposed 
into  volatile  substances  which  evaporated  under  the  laser  heating.  The 
narrow  dispersion  in  the  measured  damage  thresholds  indicates  that  this 
decomposition  occurs  at  a  well  defined  temperature. 

On  the  other  hand,  the  glass  failed  by  thermal  fracture,  as 
shown  in  Fig.  11.  It  is  possibly  because  the  surface  first  undergoes 
compressive  stress,  due  to  thermal  expansion  when  heated  by  the  laser 
beam,  then  partially  relaxes  as  a  viscous  flow  relieves  some  of  the 
stress,  at  high  temperatures,  when  the  viscosity  of  glass  drops  dra¬ 
matically  and  finally  contracts  on  cooling  to  produce  a  tensile  stress 
sufficient  to  cause  failure  by  fracture.  Temperature  increases  of 
~ 750°C  will  cause  that  type  of  failure  in  soda-lime  glass.  To 
induce  failure  in  fused  silica,  the  temperature  changes  must  exceed 
2000°C.  Such  temperatures  being  well  above  the  point  where  silica 
flows  easily,  damage  in  this  material  would  first  appear  as  surface 
distortion  due  to  melting. 

The  surface  of  glass  can  be  abraded  during  handling.  The  flaws 
thus  created  control  the  fracture  by  concentrating  stresses  and  causing 
failure  at  levels  much  lower  than  the  theoretically  predicted  loads 
calculated  from  the  forces  between  the  atoms  of  the  glass.  Therefore, 
the  distribution  of  these  flaws  is  expected  to  be  reflected  in  the 
distribution  of  failure  thresholds  for  laser  damaging  of  glass. 

Fisher  and  Hollomon  (Ref.  15)  investigated  the  theoretical  distribution 
of  stress  failure  in  glass  for  simple  and  hydrostatic  tensions  as  well 
as  unidirectional  tension  with  equal  bidirectional  compression.  They 
assumed  an  exponential  distribution  of  flaw  sizes  and  used  the  Griffith 
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FIGURE  10  -  Scanning  electron  micrograph  of  a  sample  of  PMMA  damaged 

*> 

by  a  34-us  CO^  laser  pulse  at  8.3  J/cm~.  Damage  consists 
of  many  voids  caused  by  the  evaporation  of  the  monomer 
produced  when  the  material  is  depolymeri zed  by  heating. 


I 


UNCLASSIFIED 

22 


FIGURE  11  -  Micrograph  of  the  surface  of  a  sample  of  soda-lime  glass 

~> 

damaged  by  a  34-us  CC^  laser  pulse  at  2.8  J/cnT.  Fracture 
lines  due  to  thermal  stress  are  evident.  An  above  thresh¬ 
old  irradiance  was  used  to  produce  a  large  number  of 
fractures.  Under  the  assumption  that  the  fractures  occur  at 

micro  cracks  in  the  surface,  a  minimum  density  of  15  000 
2 

cracks  per  cm  must  have  been  present  to  produce  the 
observed  results. 
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model  in  which  the  failure  strength  of  a  defect  is  inversely  propor¬ 
tional  to  the  square  root  of  its  length.  The  theory  may  be  extended 
to  include  the  stress  induced  by  laser  heating  of  the  surface.  In  this 
case,  there  is  tensile  stress  in  two  directions  parallel  to  the  surface 
but  no  stress  perpendicular  to  it.  Qualitatively,  the  results  are 
similar  to  those  obtained  in  Ref.  15  for  the  other  stress  fields. 

Figure  12  illustrates  the  frequency  distributions  for  fracture  stresses 
in  laser-heated  glass  samples  as  a  function  of  the  number  of  crack 
flaws.  There  is  a  marked  decrease  in  the  mode  and  mean  of  the  dis¬ 
tribution  as  the  number  of  cracks  increases.  More  explicitly,  the 
most  probable  failure-causing  stress  [the  stress  for  which  there  is  a 
50%  probability  of  failure)  drops  as  more  defects  are  included  within 
the  laser  heated  region.  The  dispersion  of  the  fracture  values  also 
decreases  as  the  number  of  cracks  increases.  As  Fig.  11  shows,  the 

distance  between  fractures  is  60  to  80  um,  indicating  a  crack  density 

2 

of  at  least  15  000  to  30  000  per  cm  .  For  the  12-xl2-mm  spot  sice 
of  these  experiments,  this  corresponds  to  20  000  to  36  000  defects. 

The  coefficient  of  variation  for  a  distribution  with  this  number  of 
cracks  ranges  from  0.061  to  0.065.  This  is  approximately  50%  higher 
than  the  coefficient  of  variation  we  observed  during  our  experiments 
and  it  requires  a  crack  density  of~3xl0(’  per  cm".  Since  the  distri¬ 
bution  of  flaws  is  not  the  only  factor  accounting  for  the  observed 
variation,  it  appears  that  even  more  defects  are  required  to  explain 
the  results  with  this  simple  model. 

From  measurements  of  the  fracture  strength  of  glass,  Fisher  and 

i 

Hollomon  (Ref.  15)  estimated  a  flaw  density  of ~ 1  000  per  cm"  of 
surface  area.  The  dispersion  observed  in  the  threshold  experiment 
thus  requires  a  density  of  flaws  much  greater  than  expected,  even 
when  the  15%  to  20%  standard  deviation  of  the  coefficient  of  variations 
is  considered.  During  the  glass  heating  phase,  the  compressional 


Relative  f  r equene 


UNCLASSIFIED 

24 


FIGURE  12  -  Frequency  distributions  of  fracture  stresses  for  glass 

samples  containing  1  to  10^  cracks.  The  normalized  curves 
represent  the  probability  distribution  functions  for 
failure  as  a  function  of  uniform  tensile  stress  parallel 
to  the  surface  of  the  sample.  As  the  number  of  cracks 
present  increases  the  most  probable  and  mean  failure 
stresses  decrease  and  the  dispersion  also  diminishes. 
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stresses  produce  shear  stresses  in  any  cross  section  that  is  not 
parallel  to  the  surface.  These  stresses  may  exceed  the  dynamic  shear 
strength  of  the  glass  and  produce  many  small  cracks  or  enlarge  those 
already  present,  so  that  the  distribution  of  flaws  is  significantly 
altered  by  the  time  the  tensile  stress  develops.  It  is  also  possible 
that  our  glass  samples  initially  had  more  relatively  large  crack  flaws 
since  they  were  not  protected  against  abrasions  during  handling. 

Determining  the  irradiance  level  from  the  total  beam  energy 
and  the  beam  profile  is  the  main  source  of  error.  The  absolute  cali¬ 
bration  of  the  joulemeter  is  known  to  within  ~5%  and  variations  in 
the  height  of  the  beam  profile  are~10%.  Uncertainties  in  the  latter 
could  be  reduced  considerably  by  space  filtering  the  beam  and  by 
accurately  measuring  its  profile  with  a  pyroelectric  vidicon  camera. 
These  refinements  will  be  incorporated  for  future  measurements.  We 
estimate  that  the  overall  accuracy  of  the  threshold  levels  obtained 
with  the  current  arrangement  is  15%  -  20%,  when  the  accuracy  of  the 
measuring  and  recording  instruments  and  the  calibration  of  the  energy 
monitoring  technique  are  included. 

5.0  CONCLUSION 

We  described  an  apparatus  and  a  method  for  measuring  the 
magnitude  and  the  statistical  distribution  of  the  laser  irradiations 
required  to  damage  various  glass  and  plastic  surfaces.  We  determined 
threshold  levels  for  a  50%  probability  of  damage  by ~54  us  TEA-CO.,  laser 
pulses  for  fused  silica,  common  window  glass,  clear  acrylic,  polycar¬ 
bonate,  and  polystyrene  sheet  plastics.  Significant  statistics  were 
observed  for  all  materials,  except  for  the  acrylic  plastic  whose 
dispersion  was  at  the  limit  of  the  experimental  technique.  The 
mechanism  of  laser  induced  damage  varies  for  the  different  materials. 

The  fused  silica,  polycarbonate,  and  polystyrene  specimens  failed  by 
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surface  melting,  and  have  relatively  large  dispersions,  showing  that 
their  melting  points  are  not  well  defined.  The  acrylic  plastic 
decomposed  at  a  well  defined  level  and  soda-lime  glass  failed  by 
thermal  fracture  over  a  narrow  but  finite  range.  The  number  of  surface 
flaws  required  to  explain  the  narrow  dispersion  of  threshold  irradiation 
levels  measured  in  soda-lime  glass  is  larger  than  expected  from  glass 
strength  measurements. 
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APPENDIX  A 

Maximum  Likelihood  Estimation  by  the  Probit  Method 


Under  the  assumption  that  eq.  1  represents  the  probability 
distribution  function  (p.d.f.)  of  the  damage  producing  energy  densities, 
the  probability  of  damage  occuring  at  an  energy  density,  E,  is  given  by 
eq.  3  for  the  normal  distribution  or  by 


0 

(Y-5+V/- 3)  /2\/~y 
1 


if  Y  <  5  -V~3 
if  |y-s|  <v/~T 
if  y  >  5  +v/~3 


for  the  uniform  distribution. 


In  a  batch  of  n  samples,  irradiated  at  a  level  E,  independently 
of  one  another,  the  probability  of  r  being  damaged  is 

P(r)  =  [n\/v\{n  -  r) ! ]  P  r[l  -  P]n  r  [Al] 

Consequently,  the  probability  of  a  particular  experiment  in  which  the 
number  of  damaged  sites  is  counted  for  a  variety  of  k  irradiation 
levels  to  occur  is 

P  =  n  \n.\lv.\[ni  -  1  -  P^1  **  [A2] 

i=l 

The  maximum  likelihood  method  consists  in  determining  the  values  of  the 
parameters  of  P  which  maximize  P. 

The  regression  equation  used  in  the  probit  method  was  given  in 
eq.  4  as 

y  *  5  +  {x  -  n)/cr  or 


y  =  a  +  frc 


[A3] 
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where  a  =  5  -  y/a  and  8  =  l/o. 


The  solution  to  this  equation  is  detailed  in  Appendix  II  of  Ref.  13 
and  is  found  by  solving  the  following  two  equations  for  a  and  b. 


where  a  is  the  maximum  likelihood  estimate  for  a, 
b  is  the  maximum  likelihood  estimate  for  8, 

H  -  h  '  (Pi  -  pi>’h’ 

“<  ■ 

n.  -  number  of  sample  sites  tested  at  energy  density  level  E., 
x.  =  Joulemeter  signal  volts  corresponding  to  ff., 

is  X, 

p.  -  r./n.  where  r.  is  the  number  of  damaged  sites  at  level  E., 

%  *1  Is  V  X* 

Y.  =  the  probit  corresponding  to  i,, 

“Is  "V 

-1/2  2 

and  =  (2tt)  exp[-(l\-5)  /2]  for  the  normal  distribution, 

l/2v/-3“  if  |y.  -  5 1  <V~3~ 

or  Z .  = 

t 

0  if  | Y.  -  5 |  > V  3  for  the  uniform  distribution. 

I* 


UNCLASSIFIED 

31 


Since  the  summations  in  eq.  A4  are  functions  of  Y  and  hence 
also  of  a  and  b,  the  solution  to  eq.  A4  must  be  found  iteratively. 

The  procedure  is  to  start  with  reasonable  guesses  a\,  b\  to  a  and  (3, 
calculate  Y  from  eq.  A3,  perform  the  indicated  sums,  and  solve  eq.  A4 
for  new  estimates  a 2  and  b2.  The  process  is  repeated  with  Y  =  a2  +  b2x 
and  continued  until  a  consistent  solution  is  obtained  within  the  desired 
accuracy.  This  procedure  is  performed  with  the  APL  program  PFIT  in 
which  the  criteria  used  to  test  for  convergence  to  a  solution  is 
[la.  ,  -  a.  |/a.  +  \b  .  ,  -  b  .  |  /b .]  <  10  ,  where  j  is  the  iteration 

number.  Once  a  and  b  are  found,  the  parameters  of  the  p.d.f.  are 
determined  from: 


m  =  (5  -  a)  lb 
o  =  l/b 


[A5] 


The  variances  of  the  parameters  appearing  in  the  computer  program  are 
calculated  as  follows. 

k 


=  1/  £  \w. 

Via)  = 

Vi m)  =  [V(y)  +  C m  -  X)2/S]/b2 

•KC 

=  V(b)/b 4 

V(o/m)  =  [K(y)  +  (X)2/S  J/fc4m4 


The  fiducial  limits  to  Y  are  calculated  from  Y  ±  t  V  where  t  is 

the  Student's  t  deviate  for  the  95%  level  of  probability  used.  The 
fiducial  limits  for  x  corresponding  to  a  given  probit  Y  are  given  by 
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x  *  g(x  -  x)/(l  -  g)  ±  £[(  1  -  gO^Qy)  +  (x  -  x)2 /S^]1^2 /b(l  -  g) 

where  q  =  t2/b2S 

v  xx 

Table  A1  lists  and  describes  the  important  APL  variables  found 
in  the  program.  The  listing  appears  in  Table  A2 .  Table  A3  is  a  sample 
program  output  for  the  fused  silica  case  of  Fig.  4. 
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tabu;  m 


Variables  used  in  PFIT  program 


APL  Name 


Variable  Name 


Rescript  ion 


PD 

axx 

IPG 

LL 


P 

XX 

s 


A 

JVM 

y 

p 


u 

yy 


p. 

i 


X . 

„ l 


ai>  hi 

ai* 1'  hi*l 

y 

p 


»/  n  ,w  . 

v  t  7 


p.d.f.  identifier  flag  '.V  for  Normal, 

'U'  for  uniform 
matrix  of  } r . ,  n . ,  x. { 
convergence  flag  =  0  for  no  convergence 
=  1  for  convergence 

mask  for  empty  bins  or  bins  with  no  deaths  or  no 
survivors  (i.e.  no  sites  damaged  or  all  sites 
damaged  at  a  particular  F .) 

relative  frequency  for  signal  voltage  r. 

signal  voltage 

maximum  likelihood  estimate  of  a,  the 
standard  deviation  of  the  p.d.f. 
maximum  likelihood  estimate  of  g,  the 
mean  of  the  p.d.f. 

value  of  fitted  parameters  before  iteration  i 
value  of  fitted  parameters  after  iteration  i 
estimated  probit 
estimated  probability 
p.d.f.  =  (2tt)  ~ll<-  exp. 


f -vti  -  5 r/2) 


for  the  normal  distribution  or 


*{:: 


=  1  /2\TT~  if 


s\<>rr 


otherwise  for  the  uniform  distribution 


weight  factor  where  u.  =  Z?/P.(l 

*  ill 

matrix  of  sums  in  eq.  A4  = 

r  k  k  k 

Z  n  .u . 

■  ,  i  i 
i  =  l 


L  n  .u  .x . 

..iii 

t=l 


Z  n  ,d  ,u  . 

.  .  i  is  i 

i=\ 


Z  n  ,u  .x . 

..ill 

i=\ 


k 

Z  n  .w . x 2. 
.  ,  ill 
%  =  \ 


Z  n  .v  .x  .y  , 
..it  iPiJ 
i=  1 
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APL  Name 

Variable  Name 

Description 

k  k 

XB 

* 

{Zniwixi)/tn.wi 
i=l  i=l 

k  k 

YB 

y 

(E  n.w.u  -)/Z  n.w . 
i=l  i=l 

k 

sxx 

Sxx 

Zn.w.{.x-  -  x) 2 

..til 

i=l 

k 

SXY 

Sxy 

I  »^(*.  '  x){Vi  ' 

t=  1 

k 

SYY 

syy 

£  nji.iy.  -  y)2 
i=l  *  *  1 

CHID 

- 

n.w-ip.  -  P.)2/ 2?  =  contribution  of  ith 

i  i  r  i  i  i 

point  to  chi  square  parameter. 

CHI 

X2 

chi  square  parameter  = 

.ViVPi  -  *v2/zl 

1=1 

HET 

X2/v 

heterogeneity  factor 

HU 

V 

number  of  degrees  of  freedom;  v=?c-2 

HETF 

- 

max.(l,  HET)  heterogeneity  factor  myltiplying 

variances 

VM 

V(m) 

variance  of  m 

VY 

V[Y) 

variance  of  Y 
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Description 

variance  of  o  =  oV(2>) 
k 

variance  of  a  =  1/E  n  .w  ■ 
.  ,  v  t 
t=l 

variance  of  b  =  l/Sxx 


parameter  in  Student's  t  distribution  (includes 
heterogeneity  factor)  for  calculating  confidence 
limits 


g  =  t2/b2 Sxx 


i  n 

[t!b{\-g~)][{.\-g) *  (x-x)  2/Sxx]  terms 
used  in  fiducial  limits  of  m 


x  ♦  g(x-x)/ (1-g) 
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TABLE  A2 
Program  Listings 


LI  J 

L  2  J 

[3  j 
L  4  J 


VPL  PELT  ,rJX;LL;UL  ;  S;  >' ;  X  ;  YY  ;  V ;  D ;  y ;  IF  0 ;  HE? ;  XB ;  YB \CHID ;  V 
*»  /i  /A  =  ,7  ,  L  1  J  V ,  [  .  S  3  A 

«*°v='  *’  mr-iAL  distribution ,  »y»  pot?  uhiporm  distr 

L  'S'lTIO  J  .  ' 

I  iT  V  —  3 

nTT.'l.vi’  *rr.7  4'/  Ll/ISr  O/E  DEATH,  EtiD  </ith  at  least  ode 
sd  i  vivos 


[  5  j  /■;>/pyy  5/7? 

[  o  J  L  L  *■  o  * .  i  V  a  L  1 ;  3 

C  7  j  LL  +  LL* OkhJX 

l 3 ]  p+ll/p 

L  0  J  a  -t  -I  LJ  ri  lX  [  J  ;  J 
Cl  J  J  n 

C 11 J  *SSTIHAfh  ST  ART  IDG  PARAMETERS 

[12]  fl 

[ljj  >.  S*  I  -/aX[  1  ,  p/U  ] 

[  1 4  j  93  :  V«-.  5*  +  />(AL  1 ,  P^A  ] 

[IS]  <  4  •*■  (  d  ”  •/  i  S )  ,  i  S 
[loj  .51 :  4*7.1 
C  17  J  i'—  1[lj+,1[2J*A’A 
[  1  4  J  SO  i*(PD=  ’VS*  )  /  3  'V  ,  5  J 

[10]  i  *'  PROBABILITY  MUST  SE  </  OR  U.  RETYPE.' 

[20]  PD*-  j 
[2  1]  -*50 

[22  J  Si  '.L>*Pi«)PH(  Y-S  ) 

[23]  ^(A>*l)x(p<0)x((i2)*”,5)x*".S*(r-5)*2 

[24]  -.52 

[23]  'jJ  :  P*  (  /  +  5  +  3*.5)J2x3*.S 

[26]  A2P'>1 

[27]  P*-P*P>  0 

[20]  P-( yxP<l )+P>l 

[23]  3-  (  .  5  J  *  .  5  )  x  (  3  *  .  5  )  >  |  y-  S 

[30]  52  :  ( II  /;t:a[  2  ;  J  )  x Zx^x  ( (  p/o  )  i  p)  x  (  (  p*i  )  i  l -P)  )  * .  $ 

[31]  v;.*0  SET  TO  7.-\  TO  AVOID  *0.  WEIOHT,  W,  HAS  Z  FACTOR 

n  hake  5.7 i  correct  l  ater. 


l  3  2  j 

y/-y+( 

P-P) 

i  Z  +  3 

=  0 

[  3  3  J 

A-A  +  .  X 

N  (  A  *- 

3  j  v x a x ) ,  [  i  ]v*  yy 

L  3  4  J 

-► 

(  IL'0  = 

1 )  /  5 

4 

L  3  5  j 

t 

l-.A'C 

;  3  IXXl  ;  1 

2  J 

[jo] 

-► 

(  TuL<*  /  |  " 

1  +  i  A 

A  A) /Cl 

[37  J 

/ 

F  }  -  i 

C  3  3  j 

-> 

■;i 

[30] 

A 

[43] 

.  .  ’7  ; 

T  RE 

CULT 
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L*tl  J 

A 

[4  2  J 

S’4  :  -/♦(  >-.V.l[l  1 )  *,«-3[2  ] 

2  4  3  ] 

7- .V  1  [  2  ] 

1  4  4  J 

A  b'*X[  1  ;  2  ]  iXl  1 ;  1  ] 

>5] 

y.8«-A[  1 ;  3  ]  :•  A  [  1 ;  1  J 

t4uj 

SAa^A'C  2 ; 2 ]-A£x A[ 2  ;  1 J 

[47  J 

S  A  '!*■  a  [  2  ;  3  ]  -  XU  *  A  [  1 ;  3  ] 

l'*5  J 

S/Y-(  +/5YYx5yY^YYxV)-YixA'[  1  ;3] 

[  4i  J 

C  J,  !+■*  / C'i  l D ■*■.?'/ f  x  S 7  £■*■•/ *(  P-  p)  i  0 

L  a  0  .1 

■i'-::+yi£i  ij*-  2  +  p.u 

[  5  1  J 

73  77"-  1  r  7  6T 

L  3  2  j 

/•>>9x5x(  :•  A  [  1  ;  1  ]  )  ♦  (  (  V-XB  )  *2  )  i.S’A'A 

w  D  0  ] 

A-  /  7  7F  x  (  v  A-  [  1  ;  1  ]  )  +  (  ( A  A  -  A  «  )  *  2  )  i  5  A  A' 

[54] 

ys-(s*4)  :-SAA 

[55] 

1/S  ■/*-(  (  *A'C  1 ;  1  ]  )  +  AS x  AS i- SA A  )  x  (  S:-  V)  *4 

[56] 

Y/i<-A  [  2  ;  2  ]  rA[  1 ;  1  ] 

[o7] 

VS* iSX A 

L  5  o  J 

i'f(  7£,’inF*  .  5  )  x(  (  .732  4  )  xi .  9  6  +  2 . 39  2  i.VJ-  1  . 
4.3*  V  J  =  2  )  4- 12. 7  lx  7.7=1 

032)  +  (  3.  18x,,0=  o  )  t  ( 

[  5  y  j 

S*-rxl’x  SxS  f  SAA 

[  5  0  j 

AD*V  •.'■*Sx5x.5iA[  1 ;  1  ] 

L  o  1  ] 

VS-*-(  9-SxAS)  i  l-G 

[62] 

"KIUhES’lOLDL  .50)  =  ■5,F1'J.5,";  (•3,76.5 
*  .  5  ) 

,1]).3’a(V;(7S?Fx'/.,') 

L  6  3  J 

■<■(  /Z2< 0  )  /  71 

[64] 

AZ/-(  Z’vl-Ox  v/i)*.  5 

L  0  5  J 

+  (G<l)/iJ2 

[6b] 

-(  i< X 8)  / It 3 

[67] 

HOUSE  FIDUC  [  41  LIUT  (.95):  V  >  "), 

Fd  .  5  'u '75+  7S 

[66] 

4 

[69] 

71:*  VS  FIDUCIAL  LI  UTS  AT  SKIS?  .95  Z 

o  7  S  £  ’ 

[70] 

V  4 

[71] 

,V2:,r1.95  FIDUCIAL  LI  UTS  :  T,  F7 . 5 ,  !!•«•/ 

-'3.F6.  5,  'j  (P;,F7.5 

<  ./  <  3 ,  P7  .  S , T. )!V  « ( '*9 ;  V0 ;  *3-  '40  i  A 8*  AD ) 

[72] 

“►.V  ^ 

[73] 

.73  :  *  'CUPPER  LIUT  (  .  95  )  :  V  <  '3.  .”3 . 5  *  at 

US-  -:d 

[74] 

74 :  'MAnIAHCF.  OF  TH€  ’IS  Alt  =  3.F20.19* 

<j.V 3 *  H E T F 

[  7  o  J 

'  MSI  AID  A. ID  DUVIATIOI  =  '3,F10.S.,f!  0*1, 
X  7  S )  *  .  5  ) 

FA  .  6  ,  *; )  '1  *  t  (  S ;  (  7  £ r F 

L  7  6  j 

» •»];.  ,1 TIO  sn  -  '3,F6.4,!*i  (:3,F6. 4,'3)Vm 
) 

•HU  =  rJ.Fb.3’.t(.77;S/n 

(Si  V i (I n T F * / 5  v ) * . 5 

L  /  7  ] 

[78] 

"CIETCHOS.IEITI  FACTOR  (APPLIED  IFF 

>1)  =  3,  F5 . 3  * .»  7£  " 

[7b  J 

'PARA  A  STEP  S* 

L  6  J  J 

'3  •!  =  3,  FI  0.5,  3  9  =  3,710.6*  il 

2,i  /  A 
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[81] 

'  l'i  V  A  =  .1 .  f  1  0  .  &  .  -’J 

V  6 

=  1 ,  F 1 1 . 6  • 

il  2p(  ID  T+’BSTF*  ,S)xV 

a  ,  v  a 

C  82  J 

’!';<*>  =  !i,  Fio.  6 

<Y> 

=  :'J,F10.5' 

u  Ub;  Y  3 ) 

[83] 

*  i'j  ST  =  R.F  1 0  .  o  . ;*} 

1  /ST 

=  ;rJ,  FI  2  .  10 

•a(A-Cl;l];»A-Ll;l]> 

[  84  ] 

Vj  /;*  =  I ,  F 1 0 . 5 , 7l 

SY 

=  5,Flu . 8 ' 

al  2  p  X L 1 ; 2  3] 

[35] 

,;,15*2  =  l'i,  FI  3.  S.i'J 

9  X.Y 

=  !) ,  F 1  0 . 6  ' 

ul  2  p  A  [  2  ;  2  3] 

[bo] 

'i’j.SA'X  =  i'J ,  F 1 0 . 6  ,  i!l ; 
XX; OXYiSYY) 

5  XY 

=  '1 .  F 1 0 . 6  , 

i’J;  SYY  =  !*i , F 1 0 . 6 '  it ( .9 

[87  ] 

t  f 

[38] 

t 

RXPT' • L ' 

[89]  ’  X  H  R  P  P  ROB  IT  I  <U  YY  0 

TX  i NY  A  CHI' 

[90]  " 

[91]  •  F5  .  3 , 214  ,  F5  .  2  ,  FI .  2 . 2F3 . 3  ,  PS  .  2 , 2F9 . 3  ,  F6 . 3  »  u  (  X'X  ;  L  L  /  ffV*[ 
2  ;  ]  ;££/n,i/A'C  1 ;  ] ;P;PDF  PROBIT  P;Y  ;*/xj -YY  ;«xNxXX  ;TxjxY  ;CH 
ID) 

7 


2 PRO  30  [|J]7 
V  J.*- PROB'D  X ;  B;T 

[U  A5\4£<F<y£/3:Z’F?  NOti’-IAL  DISTRIBUTION  PROBABILITY  FUNCTION 
[2]  (Pa)  ,5)p.  31933  153  “.  353553782  1.731477937  "1.821255 

9/3  1.3302741+29 

[3  J  ««-(  W*U<0  ))  +  (  *20  )xl-iM  (02  )*“.  5)  x(  +  /f*x(  *  1+  .  2316419*  |  K) 

° \ b) **- X*X i 2 
J 

OPROBITiUlV 
'tZ-PDF  PROBIT  X 
[  1  ]  -(  PDF*  •  7<yv  1  )  /.S'l ,  52 , 93 

12]  51 :  Z*-PR0BIT’l  X 

[  3  J  -0 

[  4  J  .92  :  Z- PRO  3  ITU  X 
[5]  -0 

[o  J  .53  :  l*  PR  03  ITT  X 
[  7  J  -*o 
7 


VPH03I7N  t;j]V 
'J  Y+-PR03I TU  P;T;C 
[1]  «'-i  (P>.o)-P 

[  2  ]  /’■►(  *2  x+j  )  *  .  3 

[3]  Y-T-  (  2. 5155  17  +  rx  .  8 J28S 3+  Fx  .  010328  )*(  1  +  '."*1 . 432  788  +  1’*  .  13 

32+.001303X/*) 

[4]  .  /+-5+  /x  (  "  1  )  *P<  .  3 

7 
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Vi-PPOBIX'U  P 

[1]  y-(px 2x3*.b)+D-3*.S 
V 


TABLE  A3 

Sample  Output 


FUSED  SILICA  '■ION  23  APRIL  1930 


THRESHOLDS  .  50  )  =  .  11978  (.00039) 

.95  FIDUCIAL  LIMITS  :  . 1 1 97 5 */- . 00229  ( 

VARIANCE  OF  THE  HEAR  =  .••-000  07  99097  5S8 

STANDARD  DEVIATION  =  .00917  (  .001912) 

RATIO  S/S  =  .  0 7 6 G  (  .0160) 

CHI  SQ  (  9)  s  2.489 

HETEROGENEITY  FACTOR  ( APPLIED  IFF  >1)  = 


PARA 

VETERS 

A 

=  -8.057728 

B 

= 

109.014407 

VA 

=  .014392 

VB 

r 

516.321202 

<X> 

=  .119890 

<r> 

= 

5.012060 

su 

*  105.371089 

1/5*' 

= 

.0094902692 

sx 

=  12.632999 

SI 

= 

528.126243 

SX7 

*  1.516512 

SX.I 

= 

63.528439 

SXX 

s  .001937; 

SXY 

= 

.211137;  SYY 

11745  <  V 


277 


25.506211 


.12204) 


EXPT'L 


X 

7 

R 

P 

PROHIT 

r 

HU 

YY 

HUX 

NUY 

AC 

111 

7 

1 

.14 

3.93 

4.070 

3.232 

3.94 

.350 

13.157 

.054 

114 

25 

7 

.28 

4.42 

4.397 

13.929 

4.42 

1.591 

61.250 

.005 

116 

19 

7 

.  37 

4 .66 

4.551 

11.273 

4 .67 

1.305 

51.411 

.125 

117 

16 

5 

.  '  1 

4.51 

4.724 

9.908 

4.52 

1.162 

45.805 

.418 

119 

16 

8 

.  50 

5.00 

4.898 

10.139 

5 . 00 

1 .204 

49.558 

.129 

120 

2  1 

10 

.48 

4.94 

5.051 

13.356 

4.94 

1.605 

67.466 

.165 

122 

22 

13 

.  59 

5.23 

5.215 

13.773 

5.23 

1.677 

71.821 

.003 

123 

10 

7 

.70 

5.52 

5.378 

5.042 

5.52 

.745 

32.498 

.121 

125 

20 

1  5 

.75 

5.67 

5.542 

11.435 

5.67 

1.425 

63.359 

.186 

126 

15 

12 

.80 

5.8  4 

5.705 

7.952 

5.83 

1 .004 

45.366 

.133 

128 

9 

6 

.67 

5.43 

5.869 

4.332 

5.35 

.553 
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